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1.  Introduction 


For  nanocrystalline  materials  to  be  useful,  they  must  be  sufficiently  resistant  to  grain  growth  at 
elevated  temperatures  to  retain  their  nanoscale  grain  size  (d  <  100  nm).  One  strategy  for 
stabilizing  nanocrystalline  materials  is  to  add  segregating  solutes  to  reduce  the  grain  boundary 
energy  (1-11).  Since  grain  boundary  energy  is  the  driving  force  for  grain  growth,  a  reduction  in 
grain  boundary  energy  can  impede  or  even  entirely  inhibit  grain  growth.  The  reduction  in  grain 
boundary  energy  provided  by  a  segregating  solute  is  determined  by  the  segregation  energy,  A Gseg 
(1-11).  Since  the  values  of  A Gseg  are  usually  not  available,  they  are  estimated.  Wynblatt  and 
Chatain  (12)  recently  reviewed  the  analytical  models  on  segregation  to  grain  boundaries  (GBs) 
and  surfaces  and  addressed  the  difficulty  of  meaningful  definitions  of  segregation  enthalpy, 
entropy,  and  free  energy  among  various  issues.  The  central  equation  for  all  models  is  as  follows 
for  a  binary  system: 
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GB 
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~GB 

XB 
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-exp 


A  G7e9 

RT 


(1) 


where  XgB  and  x  IIi  are  the  mole  fractions  of  component  B  (solute)  in  the  grain  boundary  and  grain 
interior,  respectively,  A Gfeg  is  the  excess  Gibbs  energy  of  segregation,  and  R  and  T  are  the  gas 
constant  and  absolute  temperature,  respectively.  The  Gibbs  energy  of  segregation  is  thus  defined 
as  follows: 


A Gseg  =  A Hseg  -  TASseg  =  (A Hseg  -  TA5“)  -  TAS^f  =  AG“  -  TAS^f,  (2) 

where  A Hseg,  A Sseg,  and  A Slsd^gal  are  the  enthalpy,  entropy  and  ideal  entropy  of  segregation, 
respectively.  Equation  1  is  thus  obtained  by  setting  equation  2  equal  to  zero,  i.e., 


AGse9  =  0,  (3) 

indicating  zero  driving  force  for  solute  segregation  to  the  interface.  With  a  given  model  for 
A Ges*y  the  compositions  in  the  grain  boundary  and  grain  interior  can  be  obtained  and  are  further 
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used  to  evaluate  the  grain  boundary  energy  of  isotropic  or  anisotropic  systems. 

Conventionally,  the  interfacial  energy  is  defined  as  the  reversible  work  needed  to  create  a  unit  area 
of  surface  (e.g.,  grain  boundary)  at  constant  temperature,  volume  (or  pressure),  and  chemical 
potentials,  i.e.,  for  an  open  system  (12).  However,  in  practical  applications,  the  interfacial  energy 
is  measured  in  a  closed  system  (i.e.,  constant  compositions).  In  the  present  work,  we  will  first 
derive  the  expression  of  grain  boundary  energy  for  a  closed  system  by  differentiating  the  internal 
and  external  variables  and  defining  the  internal  processes,  and  then  predict  the  effect  of 
segregating  solutes  on  the  grain  boundary  energy  in  binary  systems  of  44  solvents  and  52  solutes. 
This  reduction  in  grain  boundary  energy  will  be  used  to  select  nanocrystalline  alloys  in  terms  of 
thermal  stability. 


2.  Thermodynamics  of  Grain  Growth  and  Segregation 


The  combined  first  and  second  law  of  thermodynamics  of  a  system  can  be  written  as  (13) 


dG  =  —SdT  +  VdP  +  Y  dM  -  Dd £, 

i 


(4) 


where  G,  S,  V,  and  /r,  are  the  Gibbs  energy,  entropy,  volume,  and  chemical  potential  of 
component  i  of  the  system;  T,  P,  and  N,  are  temperature,  pressure,  and  moles  of  component  i 
controlled  from  surroundings;  £  and  D  represent  the  extent  of  an  internal  process  and  its  driving 
force  and  more  than  one  simultaneous  internal  processes  can  be  considered  as  shown  below.  For 
a  system  with  a  grain  boundary  area  of  A  which  reaches  a  metastable  equilibrium  under  constant 
temperature,  pressure,  and  compositions,  each  term  in  equation  4  becomes  zero,  i.e., 


dG  =  0,  (5) 

Let  us  now  consider  an  internal  process  with  the  grain  boundary  area  changed  by  dA.  Due  to  the 
composition  difference  between  the  grain  boundary  and  the  grain  interior,  there  will  be  a 
simultaneous  redistribution  of  elements,  commonly  referred  to  as  segregation.  These  two  internal 
processes  contribute  to  the  change  of  Gibbs  energy  of  the  system  as  follows  (13): 
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dG  =  ^  (v?B  ~  lli)  d'ni  +  7o dA, 


(6) 


i 


where  /if B  and  //[  are  the  chemical  potentials  of  component  i  in  the  grain  boundary  and  grain 
interior,  respectively,  dni  is  the  change  of  component  i  from  the  grain  interior  to  the  grain 
boundary,  and  y0  is  the  change  in  G  due  to  a  change  in  A  at  constant  rq,  i.e., 
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(7) 


nfB  —  ft  -  can  be  further  written  as 
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where  Ge^GB  and  G'-' ''  are  the  partial  excess  Gibbs  energy  of  component  i  in  the  grain  boundary 
and  grain  interior,  respectively.  For  a  binary  A-B  system,  we  have  driA  =  —dnB,  and  equation  6 
can  be  rewritten  as 
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where  A Gesxeg  and  A  Gseg  are  defined  as 
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It  is  evident  from  the  above  equations  that  the  grain  boundary  energy  for  a  closed  system  is 
defined  as  (14) 
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There  are  two  significant  observations  from  the  above  derivations:  (1)  AG“  is  related  to  the 
partial  quantities  of  excess  Gibbs  energy  of  grain  and  grain  boundary,  not  the  excess  Gibbs 
energies  themselves,  and  (2)  7  =  y0  when  A Gseg  =  0,  i.e.,  equation  3  and  equation  1,  meaning 
there  is  no  redistribution  of  elements  during  a  change  of  grain  boundary  area.  The  latter  is  a 
constrained  equilibrium,  because  it  poses  a  limitation  on  an  internal  process  for  a  closed  system, 
which  would  otherwise  take  place  to  reduce  the  Gibbs  energy  of  the  system.  Therefore,  the 
widely  used  equation  1 ,  commonly  referred  to  as  Langmuir-McLean  (or  Fowler-Guggenheim) 
segregation  isotherm  is  applicable  when  the  redistribution  of  elements  is  negligible  during 
measurements  of  grain  boundary  energies.  If  this  is  not  the  case,  equation  12  should  be  used  for 
grain  boundary  energy. 

There  are  many  models  on  how  to  treat  A Hseg  and  A 5“  in  equation  2  as  discussed  by  Wynblatt 
and  Chatain  (72)  including  their  dependences  on  five  degrees  of  freedom  of  grain  boundaries.  As 
the  purpose  of  this  work  is  to  compare  the  effects  of  a  wide  range  of  alloying  elements,  we  follow 
the  Wynblatt  and  Ku  model  for  surface  segregation  (72,  75)  and  assume  A S^g  =  0.  Therefore, 
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and 


a  rrelastic  =  -24tT TV AvK BG AV Ar B  ijB_  -  A Af 

3  A  p/r  a  +  4Ga'T  b 


(15) 


In  the  above  equations,  7 a,  7 b,  &a,  and  aB  are  the  free  surface  energies  and  the  grain  boundary 
areas  per  mole  of  grain  boundary  atoms  of  components  A  (solvent)  and  B  (solute),  respectively. 
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a  is  a  bond  energy  interaction  parameter  introduced  to  distinguish  the  bond  energy  in  the  lattice 
region  from  the  grain  boundary  region.  If  a  =  0,  the  expression  is  equivalent  to  the  surface 
segregation  model.  Q  is  the  regular  solution  interaction  parameter  that  can  be  approximated  as 
four  times  the  liquid  enthalpy  of  mixing  of  an  equimolar  alloy  of  the  solute  and  solvent,  A Hmix. 
A Hmix  can  be  obtained  in  the  literature  (16,  17).  Values  for  the  liquid  enthalpy  of  mixing  are 
chosen  over  those  of  the  solid  states,  because  such  values  would  include  elastic  contributions, 
which  are  considered  separately  in  the  equation,  zi  and  zv  are  the  coordination  numbers  in  and 
out  of  the  grain  boundary  plane  regions,  and  z  is  the  coordination  number  (z  —  Z[  +  2 zv).  The 
equation  for  A Hf^stlc  is  based  on  analyses  by  Friedel  (18)  and  Eshelby  (19)  with  Nav  being 
Avogadro’s  number,  KB  being  the  bulk  modulus  of  the  component  B,  Ga  the  shear  modulus  of 
the  component  A,  and  ta  and  rB  are  the  atomic  radii  of  components  A  and  B.  Equation  13 
correctly  considers  both  the  chemical  contributions  (equation  14)  to  the  segregation  enthalpy,  as 
well  as  the  elastic  enthalpy  (equation  15),  as  proposed  by  Wynblatt  and  Chatain  (12)  and 
Wynblatt  and  Ku  (15).  Recently,  Saber  et  al.  (20)  applied  a  similar  model  augmented  with  the 
Trelewicz  and  Schuh  (8)  analytical  approach  to  four  binary  systems:  Fe-Zr,  Cu-Nb,  Cu-Zr, 
Ni-W.  In  contrast  to  Saber  et  al.,  the  present  model  uses  the  concept  of  internal  state  variables  to 
define  the  metastable  state.  Additionally,  we  apply  the  present  model  to  a  much  larger  collection 
of  binary  systems  (>  1000)  using  published  experimental  and  predicted  values,  which 
necessitated  the  development  of  thermodynamic  stability  maps  to  organize  the  information  and  to 
guide  alloy  selection/development.  The  developed  maps  allow  the  visualization  of  trends 
governing  thermodynamic  stability  that  are  not  obvious  given  analysis  of  only  a  few  systems. 
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3.  Application  to  Grain  Size  Stabilization  in  Fe-Zr  and  Other  Binary 
Systems 


3.1  Grain  Boundary  Stabilization  in  the  Fe-Zr  System 

For  a  given  system  under  constant  temperature,  pressure  and  overall  composition,  the  grain 
boundary  energy  represented  by  equation  12  is  a  function  of  xBB  and  xB,  which  are  related  by 
mass  balance  based  on  the  grain  size  d  and  the  grain  boundary  model.  The  Fe-Zr  system  is  used 
to  demonstrate  the  procedure.  An  experimental  investigation  for  the  stabilization  of 
nanocrystalline  Fe  by  additions  of  Zr  solute  was  reported  in  (10).  For  bcc  Fe  with  {110}  grain 
boundary  planes,  zt  —  4  and  zv  =  2,  and  a  =  5/6.  It  is  assumed  that 
a  a  =  &b  =  cr  =  aFe  =  NAvV^  where  VFe  is  the  atomic  volume  of  Fe  (solvent). 

For  a  bilayer  grain  boundary  model  shown  schematically  in  figure  1,  one  can  make  the  following 
approximation: 


dnB  =  2  (x%B  -  x^r) 
d  A  a 


(16) 


where  T  is  the  GB  solute  excess.  With  d  being  the  volume-averaged  grain  size,  the  mass 
conservation  in  a  closed  system  with  a  mole  fraction  xq  of  Zr  added  to  Fe  is 
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where  3/d  represents  the  grain  boundary  area  per  unit  volume  for  the  spherical  grain  shape,  and 
Vm  is  the  molar  volume  for  the  solvent  Fe.  Equation  12  can  thus  be  rewritten  in  terms  of  the 
normalized  grain  boundary  energy  7/70, 
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Figure  2  shows  a  typical  plot  of  7/70  vs.  for  a  range  of  grain  size  values  at  x[}  =  0.03  and 
T  =  550  °C  with  other  parameters  given  in  table  1.  It  can  be  seen  that  each  7/70  curve  passes 
through  a  minimum  for  each  grain  size  considered.  The  equilibrium  condition  of  the  system  with 
respect  to  grain  growth  ( 6 )  is  stipulated  by  equation  5,  which  is  equivalent  to  setting  equation  18 
equal  to  zero,  i.e.,  7/70  =  0.  Figure  2  shows  that  the  curve  with  d  =  23.1  nm  (filled  red  circles) 
has  a  minimum  at  xff  =  0.26  with  7/70  =  0.  This  grain  size  value  is  designated  as  dm.  For 
systems  with  d  <  drn,  there  will  be  a  thermodynamic  driving  force  for  grain  growth  to  reduce  the 
grain  boundary  area  until  d  =  drn.  For  systems  with  d  >  dm,  there  are  two  values  where 
7/70  =  0  as  shown  in  figure  2.  However,  it  is  evident  that  if  these  two  states  were  put  together, 
atomic  diffusion  would  take  place  due  to  the  different  values  of  x%f  and  xZr  of  the  two  states 
until  both  xz^  and  xZr  become  homogeneous  in  grain  boundary  and  grain  interior,  respectively, 
and  the  grain  size  adjusts  itself  to  dm  based  on  the  model. 

Figure  3(a)  shows  the  dm  grain  size  values  as  a  function  of  temperature  for  a  range  of  x0  alloy 
contents.  At  higher  solute  concentrations  x0,  there  is  additional  solute  to  stabilize  the  grain 
boundaries.  At  the  higher  Zr  contents,  stabilization  at  a  nanoscale  grain  size  smaller  than  100  nm 
would  be  effective  up  to  temperatures  from  800  °C  to  900  °C.  There  is  an  abrupt  destabilization 
as  temperature  increases  above  a  limit,  though.  An  alternative  plot  of  the  data  in  figure  3(a)  is 
given  in  figure  3(b)  where  1  j dm  is  plotted  as  a  function  of  T  and  x0.  Complete  loss  of 
stabilization  is  revealed  as  these  curves  approach  the  dm  — >  00  limit.  Using  a  regular  solution 
model  in  the  limit  of  small  x0  values  and  fully  saturated  grain  boundaries,  Kirchheim  (6) 
suggested  that  inverse  grain  size  vs.  InT  curves  would  be  linear  with  a  negative  slope,  in 
qualitative  agreement  with  the  plot  in  figure  3(a). 

The  grain  boundary  thermodynamic  stabilization  model  used  in  the  present  work  includes  elastic 
size  misfit  energy,  AHf^gstlc,  and  chemical  bond  energy  AH^gm.  AHf^rn  is  manifested  by  the 
enthalpy  of  mixing,  A Hmix,  of  an  equimolar  Fe-Zr  liquid  phase  through  the  bond  energy 
interaction  parameter,  a;  however,  the  A  H^"n  is  dominated  by  the  magnitude  of  the  enthalpy  of 
mixing  A Hmix.  A Hfeagstlc  always  favors  grain  boundary  segregation.  A Hmix  >  0  (demixing) 
favors  grain  boundary  segregation  whereas  A Hmix  <  0  (mixing)  favors  grain  boundary 
desegregation.  The  combination  of  these  two  contributions  (A Hmix  oc  A H^m  and  A  Hfeagstic) 
dictates  the  effect  of  an  alloying  element  on  grain  boundary  energy.  The  effect  of  A Hmix  for 
Fe-Zr  alloys  is  examined  in  figure  4  for  xq  =  0.04  by  systematically  changing  the  A Hmix  values. 
For  a  hypothetical  case  of  no  chemical  effect  with  A Hmix  =  0,  drn  =  10  nm  is  predicted  at  the 
melting  point  of  Fe.  In  an  Fe-4%Zr  alloy  annealed  at  913  °C,  transmission  electron  microscopy 
and  ion  channeling  contrast  images  indicated  a  volume  average  grain  size  of  57  ±  15  nm  (9).  In 
addition,  x-ray  diffraction  data  on  the  same  alloy  as  a  function  of  annealing  temperature  are 


7 


provided  and  plotted  as  blue  stars  in  figures  3  and  4.  A  value  of  A Hmix  =  —  24  kJ  mol-1  is  in 
excellent  agreement  with  these  experimental  results  and  is  pretty  close  to  the  —25  kJ  mol-1 
obtained  from  the  literature  (see  table  1).  It  is  notable  that  the  trends  observed  for  the  Fe-Zr 
results  in  reference  (10)  lend  support  to  quantitative  predictions  for  thermodynamic  stabilization. 
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Figure  1 .  Bilayer  grain  boundary  used  for  the  regular  solu¬ 
tion  model,  zi  indicates  in-plane  bonds,  and  zv 
indicates  out-of-plane  bonds.  The  latter  can  be 
bonds  joining  atoms  across  the  grain  boundary 
(GB)  or  bonds  joining  atoms  on  the  GB  to  the 
grain  interior  (GI).  Solvent  (A)  and  solute  (B) 
atoms  are  shown  on  the  GB  bilayers. 
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Figure  2.  Normalized  grain  boundary  energy  7/70  vs.  mol  fraction  of  Zr  atoms  on  the  grain 
boundary  for  several  different  grain  sizes:  10,  15,  23.1,  30,  50  and  1  x  1012  nm.  The 
molar  Zr  fraction  xq  is  0.03  and  the  temperature  T  is  550  °C  for  this  plot.  The  plot 
intersects  7/70  =  0  at  a  grain  size  of  23.1  nm,  which  corresponds  to  the  stabilized 
grain  size,  dm. 
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(a)  (b) 

Figure  3.  (a)  Stabilized  grain  size  dm  and  (b)  inverse  grain  size  as  a  function  of  temperature  T 
for  different  Zr  molar  fractions,  xq:  0.01,  0.015,  0.02,  0.03,  0.04  and  0.05.  Experi¬ 
mental  X-ray  diffraction  (XRD)  and  transmission  electron  microscopy  (TEM)  data  for 
an  Fe-4%Zr  alloy  are  plotted  as  blue  stars  (9,  10).  AHmvx  of  —  24  kJ  mol-1  is  used  to 
show  agreement  with  the  experimental  XRD/TEM  data.  The  inverse  grain  size  shows 
a  linear  dependence  on  temperature  T  with  some  deviation  from  linearity  at  lower 
temperatures. 
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Figure  4.  Stabilized  grain  size  drn  vs.  temperature  T  for  molar  Zr  fraction  xo  =  0.04  and 
different  AHmtx  values:  0,  —20,  —24,  —25,  —26  and  —30  kJmol-1.  Experimental 
XRD  and  TEM  data  for  an  Fe— 4%Zr  alloy  arc  plotted  as  blue  stars  (9,  10).  AHmix  = 
— 24  kJmol-1  shows  excellent  agreement  with  the  experimental  XRD/TEM  data. 
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Table  1 .  Parameters  of  the  Fe-Zr  binary  system. 


Property 

Value/Units 

Equation  or  Reference 

a 

5/6 

Ct  (-intergranular  /  ( intragranular 

z 

8 

BCC,  z  —  zi  +  2 zv 

zv 

2 

BCC,  {110} 

Zl 

4 

BCC,  {110} 

7o 

0.795  Jm”2 

Ref.  (21,  22) 

7  Zr 

1.909  Jnr2 

Ref.  (21,  22) 

7  Fe 

2.417  J  m”2 

Ref.  (21,  22) 

KZr 

89.8  GPa 

Ref.  (23) 

G  Fe 

81.6  GPa 

Ref.  (23) 

vFe 

0.0118  nm3 

Ref.  (23) 

Vzr 

0.0233  nm3 

Ref.  (23) 

— 25  kJ  mol”1 

Ref.  (16,  17) 

A  TTelastic 
LJS'rLseg 

-108  kJ  mol”1 

Equation  16 

a 

31217  m2  mol"1 

NAvVpi3 

vm 

7.107  x  10~6  m3  mol”1 

NAyVFe 
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3.2  Grain  Boundary  Stabilization  in  Other  Binary  Systems 


The  same  approach  is  subsequently  applied  to  other  binary  systems  of  44  solvents  each  with  52 
solutes.  For  instance,  table  2  lists  the  solute  concentrations  required  to  thermally  stabilize  the 
binary  system  at  a  given  grain  size  (25  nm)  with  a  given  solute  for  a  particular  solvent  (in  this 
case,  Fe)  at  a  given  temperature  (0.60T^e  «  1200  K).  For  each  solute  (of  52),  the  following 
quantities  were  obtained  for  the  solvent:  the  enthalpy  of  mixing  A Hmix,  the  free  surface  energy 
7 B,  the  bulk  modulus  KB,  the  molar  volume  V^j,  and  atomic  radii  rB.  Some  of  these  are 
required  to  calculate  the  elastic  enthalpy  term  A HfeagsUc  using  equation  15.  Then,  the  content  of 
each  solute  was  changed  from  0.1%  to  10%  in  increments  of  0.1%  to  find  the  minimum  solute 
concentration  required  to  stabilize  the  system,  i.e.,  7/70  =  0  (equation  18).  Table  2  represents 
all  solutes  and  their  minimum  concentration  x0  required  to  stabilize  the  solvent  Fe  at  a  grain  size 
of  25  nm  for  a  temperature  of  0.60T^e.  The  elastic  enthalpy  A Hfeagstlc  and  enthalpy  of  mixing 
A Hmix  are  also  included  along  with  the  grain  boundary  solute  excess  T  for  the  bilayer  model, 
which  represent  reasonable  values.  This  table  does  not  include  those  solutes  that  were  not 
predicted  to  stabilize  the  system,  7/70  >  0  for  all  x0  <  0.10%.  For  stabilization  of  a  larger  grain 
size  at  the  same  temperature,  lower  solute  concentrations  are  required  (e.g.,  for  100  nm,  4.8%  Zr 
or  0.2%  Ca).  This  would  be  expected,  since  there  is  a  lower  total  grain  boundary  area  at  larger 
grain  sizes.  Similarly,  for  stabilization  of  the  same  grain  size  at  a  lower  temperature,  lower  solute 
concentrations  are  again  required.  While  there  are  a  number  of  solutes  that  would  not  be  selected 
due  to  their  high  price,  toxicity,  low  melting  temperature,  etc.,  this  model  may  shed  light  on 
prospective  binary  nanocrystalline  systems  from  a  thermodynamic  stability  standpoint. 
Additionally,  some  solutes  may  drive  the  interfacial  energy  down,  but  could  be  embrittle  the  grain 
boundary.  Hence,  this  model  is  only  one  criteria  of  several  that  would  be  important  for  selecting 
nanocrystalline  binary  systems. 

The  solutes  that  are  effective  at  stabilizing  the  nanocrystalline  grain  structure  are  functions  of 
both  the  elastic  enthalpy  and  the  enthalpy  of  mixing.  Figure  5  is  a  plot  of  the  elastic  enthalpy 
versus  the  enthalpy  of  mixing  for  both  the  stabilizing  solutes  (table  2)  as  well  as  the  solutes  that 
do  not  appear  in  table  2.  The  red  and  black  dots  denote  the  stabilizing  and  nonstabilizing  solutes, 
and  the  size  of  the  dot  for  the  stabilizing  solutes  corresponds  to  the  magnitude  of  the  minimum 
solute  concentration  required  to  stabilize  a  grain  size  of  25  nm.  Moreover,  a  convex  hull*  of  all  of 
the  systems  (dotted  line),  the  stabilizing  solutes  (red  line),  and  the  nonstabilizing  solutes  (black 
line)  is  plotted  to  delineate  the  nanocrystalline  stability  design  space.  Interestingly,  there  is  a 
noticeable  division  of  the  solutes  that  stabilize  the  nanocrystalline  grain  structure  from  those  that 
don’t.  Therefore,  figure  5  and  similar  plots  can  be  envisioned  as  ‘nanocrystalline  stability  maps.’ 

*In  computational  geometry,  a  convex  hull  for  of  a  set  X  of  points  is  the  smallest  convex  set  that  contains  X. 
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There  are  a  few  general  trends  that  emerge  from  the  nanocrystalline  stability  map  (figure  5).  The 
first  general  trend  observed  is  that  the  larger  the  magnitude  of  the  elastic  enthalpy  (i.e.,  more 
negative)  of  the  solute,  the  smaller  the  bulk  solute  concentration  required  for  stabilization. 
Additionally,  the  magnitude  and  sign  of  the  enthalpy  of  mixing  detract  from  or  reinforce 
stabilization,  with  negative  values  tending  to  reduce  stability  in  any  given  system.  Furthermore, 
stabilized  systems  that  require  higher  solute  concentrations  tend  to  have  a  lower  driving  force  for 
segregating  to  the  boundary.  Of  the  nonstabilized  systems,  there  were  a  number  of  systems  that 
reduced  the  excess  grain  boundary  free  energy,  but  did  not  reduce  it  to  zero  (recall  that  7/70  =  0 
is  required  for  stabilized  systems);  in  a  very  limited  number  of  cases,  bulk  solute  concentrations 
greater  than  10%  reduced  the  grain  boundary  energy  to  zero.  However,  in  general,  increasing  the 
solute  content  even  to  moderately  high  levels  (>5%)  is  not  always  a  good  practice  as  precipitate 
formation  is  exacerbated  at  higher  solute  contents  and  with  increased  temperature.  Experimental 
observations  generally  show  that  secondary  phases  must  be  kinetically  hindered  (24)  from 
forming,  since  they  compete  for  the  available  solute,  causing  destabilization  of  the  boundaries  and 
rapid  grain  coarsening.  Additionally,  attention  must  be  paid  to  the  extent  of  forced  solubility  that 
can  be  attained  by  a  given  processing  method,  as  this  ultimately  limits  the  amount  of  solute 
available  for  segregation  to  boundaries. 

Using  the  above  rational,  tables  for  other  solvents  that  are  often  associated  with  nanocrystalline 
materials  are  also  presented.  The  other  solvents  are  Ni,  Cu,  Al,  Mg,  Ti,  Pd,  and  W  (tables  3-9). 
Moreover,  the  same  nanocrystalline  stability  maps  are  shown  for  the  stabilizing  solutes  in  these 
systems  as  well  (figures  6-12).  A  few  general  themes  are  observed  for  these  systems: 

1.  In  a  number  of  systems,  the  majority  of  stabilizing  solutes  have  a  positive  enthalpy  of 
mixing.  For  instance,  the  majority  of  stabilizing  solutes  in  the  Mg,  Ti,  Al,  and  W  systems 
have  a  positive  enthalpy  of  mixing.  Currently  the  authors  are  unaware  of  any 
thermodynamic  reports  of  grain  growth  prevention  for  nanocrystalline  Mg,  Ti,  Al,  or  W. 
Additionally,  it  should  be  noted  that  many  of  the  predicted  stabilizing  agents  (e.g.,  Cs,  Rb, 
K,  Ba,  Na,  In,  Sn,  TI,  Sb  and  Pb),  have  been  suggested  by  Seah  to  be  grain  boundary 
embrittling  agents  (25).  Interestingly,  in  Ti,  all  of  the  stabilizing  solutes  have  a  positive 
enthalpy  of  mixing  and  there  is  a  clear  delineation  between  the  stabilizing  and 
nonstabilizing  solutes.  In  the  W  system,  a  large  variety  of  possible  stabilizing  solutes  exist 
with  the  majority  of  the  systems  requiring  a  larger  bulk  solute  content  than  many  of  the 
other  systems  examined.  Additionally,  due  to  the  high  melting  point  of  W,  high 
temperature  stability  will  suffer  from  intermetallic  formation  at  much  lower  temperatures. 

2.  A  positive  enthalpy  of  mixing  is  not  necessarily  required  for  nanocrystalline  stability.  For 
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instance,  in  the  Ni-  and  Cu-based  stabilized  systems,  an  equal  number  of  positive  and 
negative  enthalpy  of  mixing  systems  are  indicated  over  the  given  range.  Even  for  Fe,  Al, 
Mg,  W,  there  are  several  stabilizing  solutes  with  negative  enthalpies  of  mixing.  However, 
the  stabilizing  solutes  all  tend  to  have  a  greater  magnitude  of  elastic  enthalpy  in  these  cases 
to  balance  out  the  negative  enthalpy  of  mixing.  On  the  other  hand,  the  majority  of  the 
nonstabilized  systems  have  a  negative  enthalpy  of  mixing  values  and  tend  to  form 
intermetallics.  This  is  as  expected,  since  negative  enthalpy  of  mixing  values  lower  the 
segregation  enthalpy  and  the  system  orders  resulting  in  oscillations  of  the  composition;  for 
positive  enthalpy  of  mixing  values,  segregation  is  enhanced  by  forming  clusters  and/or  by 
multilayer  segregation  (26). 

3.  The  limited  experimental  data  supports  the  grouping  of  stabilizing  and  nonstabilizing 
solutes  from  the  nanocrystalline  stability  maps.  Experimentally,  there  are  few  studied 
systems  for  either  Ni  or  Cu,  for  which  the  grain  growth  was  classified  as  being  controlled 
thermodynamically.  While  a  lack  of  grain  growth  was  observed  for  nanocrystalline  Ni-W 
and  it  is  reasoned  to  be  thermodynamic  in  nature,  estimates  by  Schuh  et  al.  suggested  only  a 
60%  reduction  in  excess  grain  boundary  energy  was  possible,  and  thus  true  equilibrium  was 
not  met  (4).  Large  reductions  in  grain  boundary  energy  were  calculated  for  W  solute  in  Ni 
with  the  present  model  (^40%-60%  at  high  solute  concentration  >10%);  since  true 
thermodynamic  equilibrium  was  not  reached,  W  is  not  suggested  to  be  a  true  stabilizing 
agent  for  Ni  (table  3).  Additionally,  Cu  has  been  experimentally  added  to  nanocrystalline 
Ni  and  has  shown  minor  improvements  in  delaying  the  onset  of  grain  growth  at 
T  <  0.33 Tm  (27);  however,  in  agreement  with  our  stability  maps,  Cu  is  not  listed  as  a 
stabilizing  solute  for  annealing  temperatures  T  <  0.50 Tm. 

There  have  been  a  few  studies  which  examined  thermodynamically  stabilizing 
nanocrystalline  Cu  with  elements  such  as  Ta  and  Zr,  which  are  both  listed  as  potential 
candidates  (table  4).  Recent  atomic  simulations  on  Ta  dissolved  in  nanocrystalline  Cu  have 
shown  that  in  the  segregated  state,  Ta  can  prevent  grain  growth  even  at  high  homologous 
temperatures  >0.70 (28).  Zr  additions  to  Cu  have  shown  mixed  support  of 
thermodynamic  stability;  the  delineation  between  kinetic  and  thermodynamic  modes  of 
stability  were  unclear  (29).  The  data  in  table  4  indicates  a  high  Zr  concentration  needed  for 
stabilization.  Interestingly,  based  on  experimental  results  of  several  binary  systems 
presented  in  the  nanocrystalline  stability  maps,  most  binary  systems  that  are  effective  at 
grain  stabilization  tend  to  lie  at  the  boundary  of  the  convex  hull  for  stabilizing  solutes 
(Fe-Zr,  Cu-Nb),  and  these  are  close  to  those  solutes  that  are  not  predicted  to  be  effective 
grain  stabilizers.  On  the  other  hand,  those  systems  that  have  not  been  effective  at  grain 
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stabilization  (Cu-Zr)  tend  to  lie  on  the  interior  of  the  convex  hull. 


4.  In  some  systems,  most  solutes  tend  to  have  a  negative  enthalpy  of  mixing,  thus  limiting  the 
thermodynamic  stability  of  the  nanocrystalline  grain  structure.  For  example,  in  the  Pd 
system,  there  is  only  one  solute  (Bi)  that  can  stabilize  the  grain  size.  However,  it  is  well 
known  that  Bi  acts  as  an  embrittlement  agent  at  grain  boundaries.  Since  most  solutes  for 
Pd  have  only  negative  enthalpies  of  mixing,  this  may  help  explain  the  lack  of  stabilizing 
solutes  in  this  nanocrystalline  system. 


Table  2.  Solutes  and  their  concentrations  that  stabilize  the  grain 
structure  in  nanocrystalline  Fe  at  temperature  T  = 
0.60 T[f  =  921  °C  (1194  K)  and  grain  size  drn  =  25 
nm. 


System 

(Fe-X) 

A  TTelastic 

1  seg 

(kJmol-1) 

^HmiX 

(kJ  mol-1) 

x0 

(at%) 

r 

(/imolm~2) 

Th 

-187 

-ii 

0.6 

6.4 

Ca 

-143 

25 

0.7 

7.6 

Pb 

-135 

29 

0.7 

7.6 

Bi 

-135 

26 

0.7 

7.6 

Sr 

-123 

34 

0.8 

8.7 

Y 

-149 

-1 

0.8 

8.5 

Ba 

-120 

37 

0.8 

8.7 

La 

-138 

5 

0.8 

8.5 

K 

-50 

81 

0.9 

9.9 

Rb 

-53 

83 

0.9 

9.9 

Sn 

-123 

11 

0.9 

9.5 

Sb 

-124 

10 

0.9 

9.5 

Cs 

-46 

85 

0.9 

9.9 

T1 

-77 

31 

1.1 

11.7 

Na 

-36 

62 

1.2 

13.1 

In 

-73 

19 

1.4 

14.4 

Mg 

-52 

18 

1.9 

19.1 

Hg 

-45 

22 

1.9 

19.5 

Ag 

-31 

28 

2.0 

20.8 

Cd 

-47 

17 

2.1 

21.0 

Li 

-14 

26 

2.5 

26.0 

Sc 

-95 

-11 

2.7 

18.4 

Au 

-39 

8 

4.4 

30.9 

Hf 

-106 

-21 

5.0 

17.1 

Zr 

-109 

-25 

5.8 

16.1 
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Table  3.  Solutes  and  their  concentrations  that  stabilize  the  grain 
structure  in  nanocrystalline  Ni  at  temperature  T  = 
0.607$*  =  872  °C  (1145  K)  and  grain  size  d,m  =  25 
nm. 


System 

(Ni-X) 

A  Hel 
(kJ  mol-1) 

^Hmix 

(kJ  mol-1) 

x0 

(at%) 

r 

(/rmol  m-2) 

Th 

-202 

-39 

0.6 

6.9 

Ca 

-149 

-7 

0.7 

8.0 

Pb 

-148 

13 

0.7 

8.2 

Bi 

-145 

10 

0.7 

8.1 

Y 

-161 

-31 

0.8 

8.8 

Sn 

-137 

-4 

0.8 

9.0 

Sb 

-136 

-1 

0.8 

9.0 

Sr 

-128 

-1 

0.9 

10.1 

Ba 

-124 

0 

0.9 

10.1 

La 

-147 

-27 

1.0 

10.5 

K 

-51 

45 

1.1 

12.7 

Rb 

-54 

47 

1.1 

12.8 

Cs 

-46 

48 

1.1 

12.7 

T1 

-85 

13 

1.2 

13.2 

Na 

-38 

32 

1.5 

17.0 

In 

-82 

2 

1.7 

17.4 

Ag 

-41 

15 

2.4 

25.6 

Hg 

-51 

8 

2.5 

25.0 

Au 

-52 

7 

2.9 

28.0 

Cd 

-55 

2 

3.4 

28.1 

Mg 

-60 

-4 

6.1 

27.4 

18 


Table  4.  Solutes  and  their  concentrations  that  stabilize  the  grain 
structure  in  nanocrystalline  Cu  at  temperature  T  = 
0.607$*  =  650  °C  (923  K)  and  grain  size  dm  =  25 
nm. 


System 

(Cu-X) 

elastic 

(kJ  mol-1) 

AHmix 

(kJ  mol-1) 

x0 

(at%) 

r 

(yumol  m-2) 

Th 

-169 

-24 

0.5 

5.4 

Ca 

-139 

-13 

0.6 

6.4 

Pb 

-123 

15 

0.6 

6.5 

Bi 

-127 

15 

0.6 

6.6 

Y 

-137 

-22 

0.7 

7.3 

Sb 

-114 

7 

0.7 

7.5 

Sr 

-121 

-9 

0.8 

8.4 

Sn 

-109 

7 

0.8 

8.5 

Ba 

-118 

-9 

0.8 

8.3 

La 

-131 

-21 

0.8 

8.2 

Rb 

-53 

27 

1.1 

11.7 

T1 

-72 

15 

1.1 

11.5 

K 

-50 

25 

1.2 

12.7 

Cs 

-46 

28 

1.2 

12.7 

In 

-67 

10 

1.3 

13.2 

Na 

-36 

16 

1.9 

19.4 

Hg 

-42 

8 

2.3 

21.9 

W 

-20 

22 

2.5 

26.1 

Cd 

-42 

6 

2.7 

24.6 

Mo 

-18 

19 

2.8 

29.0 

Re 

-14 

18 

3.4 

31.6 

Hf 

-85 

-17 

4.6 

16.3 

Zr 

-91 

-23 

6.0 

14.8 

Mg 

-47 

-3 

6.8 

26.3 

U 

-62 

-7 

7.0 

22.2 

Nb 

-40 

3 

7.2 

29.0 

Ta 

-44 

2 

8.4 

28.2 

19 


Figure  5.  Nanocrystalline  Fe  stability  map  that  plots  the  elastic  enthalpy  versus  the  enthalpy 
of  mixing  for  both  the  stabilizing  solutes  (table  2,  plotted  as  red  dots)  as  well  as  the 
solutes  that  do  not  appear  in  table  2,  plotted  as  black  dots.  The  size  of  the  dot  for  the 
stabilizing  solutes  corresponds  to  the  magnitude  of  the  minimum  solute  concentration 
required  to  stabilize  a  grain  size  of  25  nm. 
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Figure  6.  Nanocrystalline  Ni  stability  map  that  plots  the  elastic  enthalpy  versus  the  enthalpy 
of  mixing  for  both  the  stabilizing  solutes  (table  3,  plotted  as  red  dots)  as  well  as  the 
solutes  that  do  not  appear  in  table  3,  plotted  as  black  dots.  The  size  of  the  dot  for  the 
stabilizing  solutes  corresponds  to  the  magnitude  of  the  minimum  solute  concentration 
required  to  stabilize  a  grain  size  of  25  nm. 
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Figure  7.  Nanocrystalline  Cu  stability  map  that  plots  the  elastic  enthalpy  versus  the  enthalpy 
of  mixing  for  both  the  stabilizing  solutes  (table  4,  plotted  as  red  dots)  as  well  as  the 
solutes  that  do  not  appear  in  table  4,  plotted  as  black  dots.  The  size  of  the  dot  for  the 
stabilizing  solutes  corresponds  to  the  magnitude  of  the  minimum  solute  concentration 
required  to  stabilize  a  grain  size  of  25  nm. 
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Figure  8.  Nanocrystalline  A1  stability  map  that  plots  the  elastic  enthalpy  versus  the  enthalpy 
of  mixing  for  both  the  stabilizing  solutes  (table  5,  plotted  as  red  dots)  as  well  as  the 
solutes  that  do  not  appear  in  table  5,  plotted  as  black  dots.  The  size  of  the  dot  for  the 
stabilizing  solutes  corresponds  to  the  magnitude  of  the  minimum  solute  concentration 
required  to  stabilize  a  grain  size  of  25  nm. 
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Figure  9.  Nanocrystalline  Mg  stability  map  that  plots  the  elastic  enthalpy  versus  the  enthalpy 
of  mixing  for  both  the  stabilizing  solutes  (table  6,  plotted  as  red  dots)  as  well  as  the 
solutes  that  do  not  appear  in  table  6,  plotted  as  black  dots.  The  size  of  the  dot  for  the 
stabilizing  solutes  corresponds  to  the  magnitude  of  the  minimum  solute  concentration 
required  to  stabilize  a  grain  size  of  25  nm. 


24 


Figure  10.  Nanocrystalline  Ti  stability  map  that  plots  the  elastic  enthalpy  versus  the  enthalpy 
of  mixing  for  both  the  stabilizing  solutes  (table  7,  plotted  as  red  dots)  as  well  as  the 
solutes  that  do  not  appear  in  table  7,  plotted  as  black  dots.  The  size  of  the  dot  for  the 
stabilizing  solutes  corresponds  to  the  magnitude  of  the  minimum  solute  concentration 
required  to  stabilize  a  grain  size  of  25  nm. 
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Figure  11.  Nanocrystalline  Pd  stability  map  that  plots  the  elastic  enthalpy  versus  the  enthalpy 
of  mixing  for  both  the  stabilizing  solutes  (table  8,  plotted  as  red  dots)  as  well  as  the 
solutes  that  do  not  appear  in  table  8,  plotted  as  black  dots.  The  size  of  the  dot  for  the 
stabilizing  solutes  corresponds  to  the  magnitude  of  the  minimum  solute  concentration 
required  to  stabilize  a  grain  size  of  25  nm. 
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Figure  12.  Nanocrystalline  W  stability  map  that  plots  the  elastic  enthalpy  versus  the  enthalpy 
of  mixing  for  both  the  stabilizing  solutes  (table  9,  plotted  as  red  dots)  as  well  as  the 
solutes  that  do  not  appear  in  table  9,  plotted  as  black  dots.  The  size  of  the  dot  for  the 
stabilizing  solutes  corresponds  to  the  magnitude  of  the  minimum  solute  concentration 
required  to  stabilize  a  grain  size  of  25  nm. 
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Table  5.  Solutes  and  their  concentrations  that  stabilize  the  grain 
structure  in  nanocrystalline  A1  at  temperature  T  = 
0.607$  =  396  °C  (669  K)  and  grain  size  dm  =  25 
nm. 


System 

(Al-X) 

A  jrelastic 
L^L1  seg 

(kJmol-1) 

^ H mix 

(kJ  mol-1) 

x0 

(at%) 

r 

(/imolm~2) 

Ca 

-97 

-20 

1.0 

7.1 

Bi 

-66 

10 

1.0 

7.6 

Rb 

-46 

25 

1.1 

8.5 

Sr 

-91 

-18 

1.1 

7.7 

Cs 

-41 

26 

1.1 

8.4 

Ba 

-93 

-20 

1.1 

7.7 

K 

-42 

23 

1.2 

9.2 

Pb 

-50 

10 

1.4 

10.3 

Sb 

-47 

2 

2.0 

13.6 

T1 

-29 

11 

2.1 

15.2 

Na 

-23 

13 

2.2 

16.1 

Sn 

-35 

4 

2.8 

18.5 

In 

-22 

7 

3.4 

23.3 

B 

-47 

0 

3.9 

20.7 
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Table  6.  Solutes  and  their  concentrations  that  stabilize  the  grain 
structure  in  nanocrystalline  Mg  at  temperature  T  = 
0.60T m9  =  389 °C  (662 K)  and  grain  size  drn  =  25 
nm. 


System 

(Mg-X) 

j^f  elastic 

(kJ  mol-1) 

^Hmix 

(kJ  mol-1) 

x0 

(at%) 

r 

(/imolm-2) 

Rb 

-38 

23 

1.2 

6.5 

Cs 

-36 

25 

1.2 

6.5 

Ba 

-64 

-4 

1.4 

7.0 

K 

-32 

20 

1.5 

8.1 

Cr 

-35 

24 

1.7 

9.2 

Sr 

-57 

-4 

1.8 

8.4 

B 

-70 

-4 

1.9 

8.8 

Mo 

-16 

36 

1.9 

10.5 

W 

-16 

38 

1.9 

10.6 

Fe 

-36 

18 

2.0 

10.6 

Nb 

-7 

32 

2.3 

12.8 

V 

-20 

23 

2.4 

13.0 

Ca 

-54 

-6 

2.5 

10.0 

Ta 

-7 

30 

2.5 

13.8 

Be 

-62 

-3 

2.6 

10.8 

Mn 

-33 

10 

2.6 

13.1 

Re 

-19 

21 

3.0 

16.3 

Na 

-11 

10 

3.6 

18.2 

Ti 

-8 

16 

3.7 

19.8 

Co 

-41 

3 

5.0 

19.7 
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Table  7.  Solutes  and  their  concentrations  that  stabilize  the  grain 
structure  in  nanocrystalline  Ti  at  temperature  T  = 
0.60 TJ!  =  996  °C  (1269  K)  and  grain  size  dm  =  25 
nm. 


System 

(Ti-X) 

A  uelastic 

LXL1  seg 

(kJmol-1) 

AHmix 

(kJ  mol-1) 

x0 

(at%) 

r 

(/rmolm”2) 

Rb 

-45 

100 

1.0 

7.4 

Sr 

-89 

53 

1.0 

7.3 

Cs 

-41 

104 

1.0 

7.4 

Ba 

-91 

57 

1.0 

7.3 

K 

-41 

94 

1.1 

8.1 

Ca 

-96 

43 

1.1 

7.9 

Na 

-22 

68 

1.5 

11.0 

La 

-73 

20 

1.8 

12.3 

Y 

-65 

15 

2.2 

14.6 

Th 

-77 

8 

2.3 

14.7 

Li 

-2 

34 

2.8 

19.8 

Mg 

-10 

16 

8.8 

21.5 

Table  8. 

Solutes  and  their  concentrations  that  stabilize  the  grain 
structure  in  nanocrystalline  Pd  at  temperature  T  = 
0.60T™  =  931  °C  (1204  K)  and  grain  size  dm  =  25 
nm. 

System 

jLf  elastic 

AH mix 

x0 

r 

(Pd-X) 

(kJ  mol”1) 

(kJ  mol-1) 

(at%) 

(ytimol  m”2) 

Bi 

-93 

-21 

6.6 

14.6 

30 


Table  9.  Solutes  and  their  concentrations  that  stabilize  the  grain 
structure  in  nanocrystalline  W  at  temperature  T  = 
0.60T^  =  2044  °C  (2317  K)  and  grain  size  dm  =  25 
nm. 


System 

(W-X) 

A  inelastic 
LX11 seg 

(kJmol-1) 

^ H mix 

(kJ  mol-1) 

x0 

(at%) 

r 

(/imolm^2) 

Rb 

-48 

129 

1.2 

9.7 

K 

-44 

124 

1.3 

10.6 

Cs 

-42 

132 

1.3 

10.6 

Ba 

-102 

74 

1.3 

10.3 

Ca 

-115 

57 

1.4 

11.0 

Sr 

-103 

70 

1.4 

11.1 

Na 

-26 

97 

1.6 

12.8 

Pb 

-87 

49 

1.8 

13.8 

Bi 

-96 

45 

1.8 

13.8 

La 

-101 

32 

2.0 

14.8 

Th 

-129 

12 

2.2 

15.2 

Y 

-102 

24 

2.3 

16.6 

T1 

-46 

52 

2.3 

17.8 

Sb 

-80 

25 

2.7 

19.3 

Sn 

-71 

27 

2.8 

20.2 

In 

-40 

38 

2.9 

21.7 

Li 

-5 

50 

3.0 

23.2 

Hg 

-22 

38 

3.3 

24.7 

Mg 

-23 

38 

3.4 

25.6 

Ag 

-2 

43 

3.6 

26.1 

Cd 

-17 

33 

4.1 

25.8 

B 

-173 

-31 

5.0 

14.3 

Sc 

-49 

19 

5.4 

25.1 

Be 

-113 

-3 

5.9 

22.1 
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3.3  Thermodynamic  Stability  Map  Changes  with  Temperature 


The  thermodynamic  stability  plots  were  also  plotted  as  a  function  of  temperature  for  solutes  in  the 
Fe  system  (figure  13).  In  general,  the  plots  do  not  change  very  much,  aside  from  a  solute  or  two. 
In  this  case,  U  is  a  stabilizing  solute  for  Fe  at  low  temperatures,  but  U  is  not  predicted  to  be  a 
stabilizing  solute  at  concentrations  below  10%  for  temperatures  T  >  0.6 Tm.  Hence,  the  bounds 
of  the  convex  hull  for  stabilizing  solutes  shift  towards  solutes  with  higher  enthalpy  of  mixing 
values  and  higher  magnitudes  of  the  elastic  enthalpy.  However,  the  bounds  do  not  change 
drastically  for  this  model. 

Additionally,  the  nanocrystalline  stability  maps  were  plotted  for  the  remaining  systems  at 
temperatures  of  0.3Tm  and  0.8Tm.  Thos  plots  are  shown  in  figures  14  and  15.  Again,  a  general 
trend  is  that  the  bounds  of  the  convex  hull  for  stabilizing  solutes  at  higher  temperatures  shift 
towards  solutes  with  higher  enthalpy  of  mixing  values  and  higher  magnitudes  of  the  elastic 
enthalpy.  Clearly,  the  bounds  for  nonstabilizing  solutes  expands  in  this  direction  as  well.  Such 
thermodynamic  stability  maps  can  be  used  for  help  in  selecting  binary  nanocrystalline  alloys  that 
are  stable  over  a  range  of  temperatures.  Materials  with  applications  at  intermediate  temperatures 
may  be  able  to  use  different  solutes  for  mitigating  grain  growth  than  applications  at  much  higher 
temperatures.  This  information  may  also  be  combined  with  density,  price,  and  embrittlement 
effect  to  aid  in  the  selection  of  solutes  for  particular  systems. 
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Elastic  enthalpy,  AH'‘”“C  (kj  mol'1)  Elastic  enthalpy,  AH*stic  (kj  mol"1) 


(c)  T  =  0.5Tm  =  1040  K 


(d)  T  =  0.6Tm  =  1194  K 


(e)  T  =  0.7Tm  =  1347 K 


(f )T  =  0.8  Tm  =  1501  K 


Figure  13.  Nanocrystalline  Fe  stability  maps  as  a  function  of  temperature:  (a)  0.3Tm,  (b)  0.4Tm,  (c) 
0.5Tm,  (d)  0.6 Tm,  (e)  0.7 Tm,  and  (f)  0.8 Tm.  As  temperature  increases,  the  amount  of  solute 
required  to  stabilize  the  grain  structure  increases,  up  to  the  point  that  some  solutes  are  effective 
at  stabilizing  at  lower  temperature,  but  are  not  effective  for  higher  temperature  stability  (e.g., 
notice  the  large  red  dot  for  U  at  low  temperatures  shift  to  a  black  dot  -  nonstabilizing  -  at  higher 
temperatures). 
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Elastic  enthalpy,  All''"1'1'  (kj  mol"1) 


(a)  T  =  0.3 Tm 


Elastic  enthalpy.  AH'‘“stic  (kj  mol"1) 


(b)  T  =  0.3 Tm 


Elastic  enthalpy,  AH^aslic  (kj  mol"1) 


(c)  T  =  0.3 Tm 


-250  -200  -150  -100  _  -50  0  50 

Elastic  enthalpy,  AH^astic  (kj  mol-1) 


(d)  T  =  0.3 Tm 


(e)  T  =  0.3 Tm 


(f )T  =  0.3  Tm 


Figure  14.  Nanocrystalline  stability  maps  at  T  =  0.3 Tm  for  a  number  of  ele¬ 
ments:  (a)  Al,  (b)  Cu,  (c)  Mg,  (d)  Ni,  (e)  Ti,  and  (f)  W. 
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Elastic  enthalpy,  Al[‘''"stic  (kj  mol"1)  Elastic  enthalpy,  Ali‘''“'Iic  (kj  mol"1) 
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Elastic  enthalpy,  AH'‘”“C  (kj  mol'1) 


(b)  T  =  0.8 Tm 


Elastic  enthalpy,  Al  l“'“stic  (kj  mol"1) 


(c)  T  =  0.8 Tm 


-250  -200  -150  -100  _  -50  0  50 

Elastic  enthalpy,  AH^astie  (kj  mol-1) 


(d)  T  =  0.8 Tm 


(e)  T  =  0.8 Tm 


(f )T  =  0.8  Tm 


Figure  15.  Nanocrystalline  stability  maps  at  T  =  0.8 Tm  for  a  number  of  ele¬ 
ments:  (a)  Al,  (b)  Cu,  (c)  Mg,  (d)  Ni,  (e)  Ti,  and  (f)  W. 
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3.4  Relationship  between  Solute  Concentrations 


The  relationship  between  the  bulk  solute  concentration,  the  grain  boundary  solute  concentration 
and  the  grain  interior  solute  concentration  for  all  stabilized  solvent-solute  combinations  (566 
stabilized  systems  total  from  1294  total  systems)  is  discussed  herein.  There  are  a  few  trends  that 
appear  from  examining  these  relationships.^  As  the  bulk  solute  concentration  required  to 
stabliize  the  nanocrystalline  grain  structure  increases,  the  concentration  in  the  grain  boundary 
increases  as  well  (figure  16).  In  figure  16(a),  notice  that  the  grain  boundary  concentration 
increases  dramatically  as  a  function  of  this  minimum  bulk  solute  concentration;  this  trend 
suggests  that  there  is  a  limit  to  the  amount  of  bulk  solute  concentration  that  can  be  added  to 
stabilize  some  systems.  In  figure  16(b),  the  grain  interior  solute  concentration  deviates  from  the 
1:1  line  drawn  as  solute  segregates  from  the  grain  interior  to  the  grain  boundary,  effectively 
resulting  in  a  lower  grain  interior  concentration.  As  the  stabilized  grain  size  increases,  and  grain 
boundary  volume  decreases,  the  solute  concentration  in  the  grain  interior  will  approach  the  bulk 
solute  concentration,  i.e.,  less  bulk  solute  concentration  is  required  to  stabilize  the  grain  size. 

The  ratio  between  the  solute  concentration  in  the  grain  boundary  and  the  grain  interior  changes  as 
a  function  of  the  bulk  solute  concentration.  Figure  17(a)  plots  the  bulk  solute  concentration 
against  the  ratio  of  solute  in  the  grain  boundary  to  the  grain  interior  for  the  stabilized  systems  at  a 
grain  size  of  25  nm  and  a  temperature  of  0.6 Tfy  (T^r  is  the  melting  temperature  of  the 
corresponding  solvent  A).  As  the  bulk  solute  concentration  to  stabilize  the  nanocrystalline  grain 
structure  increases,  the  ratio  of  solute  concentration  in  the  grain  boundary  to  the  grain  interior 
decreases  (figure  17).  The  solute  concentration  ratio  axis  in  figure  17  is  plotted  on  a  log  scale  to 
highlight  the  large  difference  in  this  ratio  as  a  function  of  bulk  solute  concentration.  This  sort  of 
trend  is  as  expected,  i.e.,  the  systems  that  stabilize  at  lower  bulk  solute  concentrations  do  so 
because  of  a  greater  driving  force  for  solute  to  segregate  to  the  grain  boundary  and,  hence, 
possess  a  higher  ratio  of  solute  in  the  boundary  than  the  grain  interior.  Also,  notice  that  the  solute 
concentration  ratio  in  this  model  is  >10-100  times  higher  in  the  grain  boundary  than  in  the  grain 
interior  for  most  systems,  which  may  dramatically  impact  properties.  Figure  17(b)  shows  these 
same  systems  for  several  different  stabilized  grain  sizes.  First,  as  the  stabilized  grain  size 
increases,  the  ratio  of  solute  concentration  in  the  grain  boundary  to  the  grain  interior  increases 
dramatically  for  those  systems  that  require  a  low  bulk  solute  concentration  for  stabilization.  This 
result  has  to  do  with  the  rapid  decrease  in  grain  boundary  volume  with  increasing  grain  size 
combined  with  a  reduced  amount  of  bulk  solute  required  to  stabilize  the  grain  structures,  which  is 

'  There  are  a  few  systems  that  appear  to  be  outliers  from  the  general  trends.  For  instance,  the  two  systems  above 
the  1:1  line  in  figure  16(b)  are  the  K-Pt  and  Cs-Pt  systems.  In  figure  17,  there  were  multiple  systems  that  were 
removed  as  outliers:  K-Pt,  Cs-W,  Cs-Re,  Cs-Os,  Cs-Ir,  Cs-Pt. 
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more  pronounced  in  the  dilute  limit.  By  comparison,  at  larger  bulk  solute  concentrations,  there  is 
a  much  smaller  difference  in  this  ratio  as  a  function  of  the  stabilized  grain  size. 
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(a) 


(b) 

Figure  16.  (a)  Grain  boundary  solute  concentration  XgB 
and  (b)  grain  interior  solute  concentration  xJB 
as  a  function  of  the  minimum  bulk  solute  con¬ 
centration  required  to  stabilize  the  nanocrys¬ 
talline  grain  size  at  a  grain  size  of  25  nm  and  a 
temperature  of  0.6T('f  (TBf  is  the  melting  tem¬ 
perature  of  the  corresponding  solvent  A). 
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Figure  17.  The  ratio  of  the  solute  concentration  in  the 
grain  boundary  to  the  grain  interior  as  a  func¬ 
tion  of  the  bulk  solute  concentration  for  stabi¬ 
lized  nanocrystalline  systems  at  a  temperature 
of  0.6 Tfy  for  (a)  a  stabilized  grain  size  of  25  nm 
and  (b)  stabilized  grain  sizes  of  10,  25,  50  and 
100  nm. 
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4.  Summary  and  Conclusions 


A  fundamental  thermodynamic  analysis  is  performed  for  the  change  in  grain  boundary  energy 
that  results  from  solute  segregation  to  grain  boundaries  in  a  closed  system.  When  the  grain 
boundary  energy  is  reduced  to  zero,  the  grain  size  can  be  stabilized  against  further  change.  The 
theoretical  results  provide  a  detailed  picture  of  the  stabilization  effects.  The  developed  model 
and  stability  maps  aid  in  selecting  binary  systems  that  will  open  new  regimes  in  processing, 
manufacturing  and  applications  space.  Specific  findings  of  the  present  work  include: 

1.  A  methodology  and  model  for  thermodynamic  stabilization  of  nanocrystallinity  is 
presented.  This  thermodynamic-based  model  allows  for  the  prediction  of  a  large  number 
of  individual  binary  systems. 

2.  Numerical  calculations  were  made  for  the  Fe-Zr  system,  demonstrating  that 
thermodynamic  stabilization  for  nanoscale  grain  sizes  can  be  expected  in  Fe-Zr  alloys  up  to 
900  °C  for  Zr  contents  on  the  order  of  4  at%  Zr  (figures  3-4).  This  is  in  good  agreement 
with  the  experimental  observations. 

3.  A  thermodynamic  stability  map  was  made  for  Fe-based  systems.  A  strong  effect  of 
enthalpy  of  mixing  is  revealed  in  the  calculations  (figure  5).  We  conclude  that  in 
appropriate  alloy  systems  with  suitable  values  of  elastic  and  mixing  enthalpy  of 
segregation,  thermodynamic  stabilization  of  a  nanoscale  grain  size  is  possible  up  to  high 
homologous  temperatures. 

4.  The  utility  and  universality  of  this  model  are  demonstrated  by  applying  it  to  over  2000 
binary  systems.  Herein,  we  applied  this  to  Fe,  Ni,  Cu,  Al,  Mg,  Ti,  Pd,  and  W  (tables  2-9) 
However,  these  maps  must  be  balanced  with  other  properties,  e.g.,  secondary  phase 
formation,  solubility  limit,  and  grain  boundary  embrittlement  must  be  considered  as  well. 
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